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ABSTRACT 

We report the discovery of a hot Jupiter transiting a subgiant star with an orbital period of 6.87 days thanks to public photometric 
data from the Kepler space mission and new radial velocity observations obtained by the SOPHIE spectrograph. The planet KOI-428b 
with a radius of 1.17 ± 0.04 Rj„p and a mass of 2.2 ± 0.4 Mj„p , orbits around a F5IV star with i?* = 2.13 ± 0.06 Rq , M*= 1.48 ±0.06 
Mq and Te(f = 6510 ± 100 K. The star KOI-428 is the largest and the most evolved star discovered so far with a transiting planet. 

Key words, planetary systems - techniques: photometry - techniques: radial velocities - techniques: spectroscopic, star : 
individual(KOI-428, KIC10418224, 2MASS 19471528+4731357) 



1. Introduction 



2. Kepler observations 



Launched in march 2009, Kepler is the second space mission 
designed to find transiting exoplanets with high accuracy pho- 
tometry. This mission has akeady shown its capability to dis- 
cover new planets during the first 33.5 days of science operations 
with the announcement of a hot-Neptune-hke planet (Borucki et 
al. 2010a^, four hot-Ju piters (Koch et al. |2010; Je nkins et al., 



2010a| [Dunham et al.| [20T0 



Latham et a l. 2010), and a sys 
tem of two transiting Saturns ( Hoknan et al. ||201()| . From the 
156,000 stars observed, 706 exoplanet candidates were identi 
fied and 306 of them were published by Borucki et al.|(|2010bll. 



All these candidates are detected around faint stars (14 < mV 
< 16) not included in the ground-based follow-up conducted by 
the Kepler team. 

To establish the planetary nature of these candidates and to 
assess the fraction of false positives, a high-resolution spectro- 
scopic follow-up must be carried out. This could allow various 
configurations of eclipsing binaries to be discarded and the prop- 
erties (mass, density) of actual planets to be established ( |Bouchy| 
|et al.|[2009b[ ). We present here the radial velocity follow-up per- 
formed with the SOPHIE spectrograph, mounted on the 1.93- 
m telescope at Observatoire de Haute Provence, of the Kepler 
Object of Interest KOI-428 as part of a first set of Kepler tar- 
gets followed-up with SOPHIE (Bouchy et al., in prep). This led 
to the discovery and the characterization of the new exoplanet 
KOI-428b. 



Send offprint requests to: A. Santeme 
e-mail: alexandre . santerne@oainp . fr 

* Based on observations made with the 1.93-m telescope at 
Observatoire de Haute-Provence (CNRS), France 



The target KOI-428 was observed by Kepler duiing the first 
33.5-day segment of the science operations (Ql) from May, 13 
to June, 15, 2009 with a temporal sampling of 29.4 min. The 
various identifications of this target, including coordinates and 
magnitudes are listed in Table [l] The publicly-avaiIabl^A'£7?/er 
light curve contains 1639 points, thirteen of which are affected 
by known systematic effects, such as the loss of fine pointing 
caused by momentum desaturation or the occurrence of strong 
argabrightenin g^ and were discarded by the Kepler pipeline 
Penkins e t al. 2010b). Based on the four transits that occurred 
during Ql, Borucki et al. (2010b) report that this system consists 
of a ~ 1 .93 solar radius star with a 1 .04 Jupiter radius companion 
in a 6.87-days orbit. Figure [T| presents the corrected light curve 
from the MAST database showing four periodic transit events 
with a depth of about 3.9 mmag. One can see from this figure 
that the first event is slightly deeper than the other three. We 
therefore studied if other known systematic effects were present 
during the transits, in particular during the first one. We found 
that none of the transits were affected by Reaction Wheel Zero 
Crossings or by a fainter argabrightenings, which were not dis- 
carded by the pipeline. Also, no difference was observed in the 
out-of-transit flux of the first transit as compared to the other 
three, which would have indicated that the difference in depth 
of the first transit was caused by the detrending performed by 



' http : //ar chive . stsci . edu/kepler/data_search/ 

Jsearch.php 

- Name given to cadences for which all the focal plane is illumi- 
nated by a presently unexplained effect. For more informations, see 
http : //archive . stsci . edu/kepler/release_notes/re lease_| 
|notes5/Data_Release_65_2(S18(S6Q414.pdf, Section 6.1 
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Table 1. KOI428 IDs, coordinates and magnitudes 



Kepler Input Catalog (KIC) 


10418224 


Kepler Object of Interest (KOI) 


428.01 


2MASS ID 


19471528+4731357 


Right Ascension (J2000) 


19:47:15.29 


Declinaison (J2000) 


+47-31-35 8 


Kepler magnitude'"' 


14.58 ± 0.02 




14.9 




14.3 




13.536 ± 0.021 




13.294 ± 0.022 




13.27 ±0.30 


E{B - 


0.139 ± 0.1 



from the Kepler Input Catalog. from USNO- 
A2 catalog. from 2MASS catalog. 
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Fig. 1. Light curve available on MAST database showing the 
four transits detected on KOI-428. The red curve shows the best 
fit of the transit. 



the pipeline. On the other hand, it is known that during QO and 
Ql a couple of variable stars -an eclipsing binary and an intrin- 
sic variable- were chosen as guide stars, and that this affected 
the spacecraft guiding. As a consequence, the time-series of the 
column position of KOI428 presents positive sharp transit-like 
features lasting about 10-15 cadences. One of these "bumps" co- 
incides with the first transit and might explain the amplitude dif- 
ference observed. Since the inclusion of this transit in our fitting 
procedure does not change our results significantly, we decided 
to use the whole Kepler light curve as is provided in the MAST 
database. 



3. SOPHIE observations 

We performed the high-resolution spectroscopy follow-up of 
KOI-428 as part of a first set of Kepler candidates (Bouchy et al. 



in prep.) with the SOPHIE spectrograph ( |Pemichot et al. 2008 
Bouchy et al. 2009a[ l installed on the 1.93m telescope at the 
Observatoire de Haute-Provence, France. SOPHIE is a cross- 
dispersed high-resolution fiber-fed echelle spectrograph, stabi- 
lized in pressure and temperature and calibrated with a Thorium- 
Argon lamp, mainly dedicated to measure precise radial veloc- 
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Fig. 2. Phase-folded radial velocity curve with our best-fit (top 
panel) and residuals (bottom panel). Open circles represent the 
four observations strongly affected by the scattered Moon light. 



ity (RV) on solar-type stars for exoplanets and asteroseismology 
studies. 

Observations were made with the high efficiency mode, with 
a spectral resolution of 39,000 at 550 nm, and the slow read-out 
CCD mode. Our observations made use of both of the SOPHIE 
fibers : fiber A was centered on the target, while fiber B acquired 
the sky background ~ 2 arcmin from the target. The typical in- 
trinsic stability of SOPHIE (3 ms^'h"') does not require to use 
the simultaneous calibration and the sky fiber is critical to re- 
move the scattered Moon light. Eleven spectra of KOI-428 were 
obtained from July, 15 to September, 13, 2010 with an expo- 
sure time of 1 hour. SOPHIE spectra were reduced with the on- 
line pipeline. Radial velocities were obtained by computing the 
weighted cross-coiTelation function (CCF) of the spectra using a 



numerical spectral mask corresponding to a G2V star ( Baranne 



|et al. 1996 Pepe et al. 2002 1. Four spectra were strongly af- 



fected by the Moon background light and were corrected with 



the same procedure as for previous planets ( PoIIacco et al. 2008 



Barge et al. 200^ 



2008 1, and as described for 



jHebrard et ak] 

HARPS ( |Bonomo et al.||2010| l. This effect cannot be removed 
completely on such a faint star, hence we conservatively dou- 
bled the uncertainties of the four coiTesponding measurements 
to account for possibly remaining systematics and to get a re- 
duced close to one. Radial velocities are listed in Table|3]and 
plotted in Fig. |2] They show a clear variation in phase with the 
Kepler ephemeris listed in Table|2]and compatible with the reflex 
motion of the parent star due to a planetary companion. Since no 
significant eccentricity is seen in the data (e = 0.05+0.21) we de- 
cide to fit a Keplerian circular orbit using the updated ephemeris 
listed in Table |2] This does not affect the radial velocity semi- 
amplitude within 1-cr level. The best-fit has a semi-amplitude K 
— 179 + 27 ms"' and a cto-c - 65 ms"' which is comparable 
to the mean radial velocity uncertainty. 

To assess the possibility that the RV variations are due to a 
blended binary, we checked both the variations of the bisector 
span of the CCF as well as the dependences of the RV variations 
with different cross-correlation masks (Bouchy et al. 2009b[ l. 
The CCFs were computed with a FOV, G2V and K5V masks for 
all the spectra and no significant differences were found in the 
amplitude of the RV variations (Kfoy - 177 ± 44 m s"', Kq2v - 
\19 + 21v(\s-\Kk5v = 194±32ms"'). The bisector span, listed 
in Table [3] and plotted in Fig. [3] does not reveal any significant 
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Fig. 3. Bisector variations as function of the orbital phase (left 
panel) and as function of radial velocity (right panel). Open cir- 
cles display the four observations affected by the background 
Moon light. Bisector error bars are taken as twice the radial ve- 
locity uncertainties. No significant slope is visible between ra- 
dial velocity and bisector, discarding a blend scenario. 



variations at a level more than two times smaller than the RV 
changes. Moreover, neither radial velocities, nor their residuals 
present any significant correlation with the bisector (We find a 
Spearman's rank-order correlation coefficient between RV and 
bissector of -0.17 + 0.31). These tests allow us to secure the 
planetary nature of KOI-428b. 



4. Spectral classification 

We performed the spectroscopic analysis of the parent star 
KOI-428 using SOPHIE spectra. Individual spectra have a low 
signal-to-noise, which does not allow a proper spectral analysis. 
Therefore, we co-added order per order the first 5 spectra not 
affected by the Moon light once set to the barycentric rest 
frame. We finally get a co-added spectrum with a S/N ^ 82 
per element of resolution on the continuum at 5550 A with a 
spectral resolution of 39,000. We derived the star's photospheric 
parameters using the semi-automatic software package VWA 
( |Bruntt et aL\ [2^ [2008] ) as described in Bruntt et al.| ( |20T0| l. 
We derived vsini- 9.0 + 2 kms Vmicm = 1.0+1 kms~ , 
Teff = 6510 ± 100 K and [Fe/H] = O.rjf ■ We estimated the 
surface gravity by fitting synthetic spectra to the wings of 
the Mgl b and Nal D lines adopting the abundances derived 
from weak lines. We found log g= 4.1 + 0.2, a value that the 
quite low S/N combined to the moderate spectral resolution 
did not allow us to accurately estimate. The derived stellar 
parameters are reported in Table |2] Our stellar parameters 
are significantly different from those given by Borucki et al. 
( |2010b| ) (TefF =6127K and log §=4.55). However these authors 
note that "spectroscopic observations have not been made for 



the released stars, so it is important to recognize that some of 
the characteristics listed for the stars are uncertain, especially 
surface gravity." Using eq. 3 in Casagrande et al. (2010) with the 
J and K magnitudes from the 2MASS archive (listed in Table [T 



and a reddening E(J-K)=0.074 (computed from Cardelli et al 
1 1989) we have computed a photometric temperature TgfF 



: 6604 

+ 132 K, which is in good agreement with the spectroscopic 
value. 



5. System parameters 

The Kepler light curve was modeled using the analytical expres- 
sions from'Mandel & Agol ( 2002) with a linear limb-darkening 
law. We adjusted the radius ratio k = Rp/Ri,, the impact param- 
eter b, the system scale a//?*, and the ephemerides parameters 
P and Tfi-. The eccentricity of the orbit was set to zero in the 
light-curve analysis, since there is no compelling evidence from 
the radial velocity measurements that the orbit is eccentric. The 
coefficient of the limb darkening law was not adjusted, but was 
allowed to vary within the limits defined by the spectral analy- 
sis of Sect. |4] (see below). On the other hand, the coarse sam- 
pling of the Kepler Long Cadence light curves has been shown 
to produce an overestimation of the impact parameter ( pKipping| 

& Bakos , 2010) and hence a smaller measured stellar density. To 
1 1 i I 1 J — 1 — . 

avoid this, we follow the solution proposed by Kipping ( 2010| l, 
i.e. we generate an over sampled model light curve which we bin 
to match the sampling rate of the Kepler data before comparing 
it to the actual measurements. For the present analysis our model 
light curve was five times as densely sampled as the Kepler data. 

The best-fit parameters were found by ;^'^-minimization us- 
ing the downhill simplex algorithm (Nelder & Mead 1965), as 
implemented in the SciPy librarj]^ The fit was repeated starting 
from 1500 randomly-generated starting points, drawn from uni- 
form distributions that allowed a wide range of values for each 
parameter Additionally, the value of the limb-darkening coef- 
ficient was drawn at each iteration from a normal distribution 
centered on the value obtained by linear interpolation of the ta- 
ble by Sing (2010) using the stellar parameters computed with 
VWA (see Sect. [4]i, and standard deviation equal to half the dif- 
ference between the maximum and minimum values allowed by 
the 1-cr intervals of the stellar parameters. In this way, we try to 
take into account the systematic effects introduced by fixing this 
value during the fit. The best-fit model is shown in Fig[T]and 4] 
and the parameters are listed in Table [2] We also repeated the 
fit leaving out the first transit, which is slightly deeper than the 
other three (see Sect. [2]). We found that k is within the 1.2-cr, and 
the other parameters are all within 1-cr of the values determined 
using the entire dataset. 

To estimate the uncertainties in the obtained parameters, we 
constructed synthetic data sets scrambling the residuals of the 
best-fit model and adding them back into the model light curve. 
We repeated the processes 2000 times, refitting the data for each 
realization, and recording the obtained parameters. Here again, 
the starting points were randomly varied and the limb-darkening 
coefficient was drawn from the same distribution. If we assume 
that the set of transit parameters obtained this way is a good 
approximation of the true distribution of the best-fit parameters 
obtained with the downhill simplex algorithm, then we can use 
it to estimate the confidence intervals of these best-fit parame- 
ters. The error bars reported in Table |2] correspond to the 68% 
confidence interval, defined so that the cumulative probability 
is 16% above and below the upper and lower confidence limit, 
respectively. 

Finally, the stellar density obtained from the transit fit was 
used to calculate the value of the surface gravity (see below), 
which provided a new value for u from which we computed 
a new value for the stellar density by repeating the whole fit- 
ting process. This procedure was iterated until convergence was 
reached, for u - 0.548!|5|5|o and the values reported in Table|2] 

We used the photospheric parameters from the spectral anal- 
ysis and the stellar density from the transit modeling to deter- 



See the SciPy Web site at http: //www. scipy.org 
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Fig. 4. Unbinned phase-folded light curve and residuals of the 
four transits of KOI-428 showing the best fit model of the transit. 



mine the star's fundamental parameters in the (Tgff , m]/^ /Ri,) 
space. The location of the star in the H-R diagram was com- 
pared to STAREVO L evolution tracks (Palacios, private com.: 



Turck-Chieze et al^[201 0) minimizing the usual statistics, but 



taking into account the time, t [Gy], a given model spends in the 
error box: 



2 

^ model 



1 /Teff-T, 



eff.model 



1/3, 



model 



1 / [Fe/H] 



[Fe/H] model 



O" [Fe/H] 

We found two distincts sets of solutions which correspond to 
different evolutionary states: either the star is in the last main- 
sequence evolution stage or it is in the very beginning of the 
post main-sequence phase. The main-sequence solutions corre- 
spond to a star with mass of M* = 1.62 + 0.10 M© which is 
close to the red hook at main-sequence turnoff. A similar situa- 



tion is reported for another planet host-star, Kepler-4 ( Borucki et 
al. |2010b I whose location in the H-R diagram is also compatible 
with the two evolutionary phases that could not be distinguished. 
In the case of KOI-428, the post main sequence solutions are 
clearly separated from main sequence solutions and correspond 
to a less-massive older star M* = 1.48 ± 0.06 M© (see Figure |5]l. 
These solutions appear as the most probable due to the rapid 
evolution of stars at the end of the main sequence. We therefore 
adopted the post main-sequence solutions. While the mass for 
the two evolutionary phases differ by 14%, both sets of solutions 
give a radius that differs by less than 4% and thus have nearly no 
impact on the size of the planet. The inferred surface gravity is 
log g - 3.94 ± 0.32, in agreement with the spectroscopic value 
at 1-cr. Clearly a better estimate of the metallicity with a much 
higher quality spectrum would allow to clear out this uncertainty 
on the evolutionary status of the star and further on the mass. The 
stellar properties corresponds to a spectral type F5IV. With the 
adopted stellar parameters, we derived for the transiting planet 
Up =2.2 + 0.4 Mjup and R,, = 1.17 + 0.04 Rjup . 

6. Discussion and conclusion 

The hot Jupiter KOI-428b is the transiting planet with the largest 
host star detected so far (2.13 ± 0.06 R© ). Figure |6] shows the ra- 




M = 1.62,U , [Fe/HI = 0.1 

M = 1.48A/ ., [Fe/H] = 0.1 

M = 1.4A/,,, [Fe/H] = 0.0 

M = 1.7.1J , [Fe/Hl=0.25 



Fig. 5. STAREVOL evolutionary tracks models in the (Tgff , 
IR-i,) space for main-sequence stars (solid lines) and post- 
main sequence stars (dashed-dot lines). The purple point repre- 
sents the observed position of the star KOI-428 with the associ- 
ated error bars in this space. The two best models are the green 
solid line with = 1 .62M0 and the adopted one in dashed-dot 
red Hne with M^, - 1.48Mo (see text). 



dius of transiting planets as a function of the radius of their host 
stars. Only 5 transiting planets orbit a star with a radius larger 
than 1.8 Rq including Kepler-5 (Koch et al. 2010|l, Ke pler- 
7 ( [Latham et aTl[2 0T0'), TrES-4 (Ma ndushev et al.| |2007| ) and 
HAT-P-7 (|Pal et al. , 2008). Transiting planets around such large 
stars are difficult to detect from ground-based photometry due 
to the long transit duration and the small transit depth. We note 
that only 8 candidates over the 306 published in Borucki et al. 
( |2010 b) have an estimated host star radius larger than 1.8 Rq , 
including KOI-428. 



From our determination of the stellar radius and the pro- 
jected rotational speed of the star i/sin/, we note that the mini- 
mum star rotational period is Prot„„„ - 11.9 + 2.7 days, which 
is compatible with a synchronization at twice the planet or- 
bital period. This is also the case of other transiting exoplan- 
ets (e.g. X O-4b ( |McCullough"era l. 2008), HAT-P-6b ( [Noyes et 
[aL] [2058) 1, HAT -P-8b ([Latham et al.i|2009) ) with host-star hot- 
ter than 6000K ( [Lanza 2010| l. Additional transit measurements 
from Kepler will permit to determine the true stellar rotational 
period (as for CoRoT-4b (fAigrain et al. 



2008 Moutou et al. 



|2008[ )) and permit to understand the star-planet interaction anc 
evolution. 

Planetary mass error bars are dominated by radial ve- 
locity uncertainties and stellar characterization uncertainties. 
Additional higher S/N spectra which permit a better estimation 
of [Fe/H] and Tetf will definitely reduce the star mass uncertain- 
ties and the discrepancy between evolution solutions. 

This detection demonstrates the efficiency of small tele- 
scopes with dedicated instrumentation for the ground-based 
follow-up of space missions like Kepler and CoRoT ( [Baglin et 
|al. 2006 ). The equivalent of 1 night on a 2-m telescope led to 
estabUsh the planetary nature of a transiting candidate and to de- 
termine its mass. Only the follow-up of a significant amount of 
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Table 2. Star and planet parameters. 



1%. 



a - 



Ephemeris 


Planet orbital period P [days] 


6.87349 ± 0.00064 


Transit epoch T„ [BJD] 


2455005.5198 ±0.0024 


Results from radial velocity observations 




llrnitnl *=>r'r'f»nti"ir'it\/ *3 




Semi-amplitude A" [ m s"'] 


179 ± 27 


Systemic velocity V, [kms"'] 


6.565 ± 0.020 


O-C residuals [ms"'] 


65 


Fitted transit parameters 




Radius ratio k = Rp/R^ 


0.0565 ± 0.0004 


Scaled semi-major axis a/R^, 


8.11 ±0.06 


Impact parameter b 


0.045!°;S3J 


Deduced transit parameters 




M^'V/^* [solar units] 


0.533 ± 0.004 


Stellar density [^ c'i?7 "^] 


213 + 005 


Inclination / [deg] 


89.7!°^ 


Transit duration Tin [h] 


6.86 ± 0.06 


Spectroscopic parameters 




Effective temperature Teff [K] 


6510 ± 100 


Metalhcity [Fe/H] [dex] 


10+°" 


Stellar rotational velocity u sin ; [ km s"' ] 


9.0 ±2 


Spectral type 


F5IV 


Stellar and planetary physical parameters from combined analysis 


Star mass [Mq ] 


1.48 ±0.06 


Star radius R^, [Rq ] 


2.13 ±0.06 


Surface gravity log g" [dex] 


3.94 ± 0.32 


Age of the star t [Gy] 


2.8 ±0.3 


Distance of the system [pc] 


2700 ± 200 


Orbital semi-major axis a* [AU] 


0.080 ± 0.003 


Planet mass M,, [M j ] 


2.2 ±0.4 


Planet radius i?,,[Rj] 


1.17 ±0.04 


Planet density pp [g cm^^] 


1.68!-; 


Equilibrium temperature ^ T^,, [K] 


1620 ± 30 



" derived from and R* . 
* derived from a/R* and R* . 
considering a zero albedo and a perfect atmospheric thermal circu- 
lation. 



Kepler candidates will provide the real fraction of true transiting 
exoplanets. 
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Fig. 6. Radius-radius diagram of all transiting exoplanets discov- 
ered so far, overplotted with curves of equal transit depth. Open 
circles are transiting planets discovered by ground-based pho- 
tometry, filled circles are those discovered by space-based pho- 
tometry and the red square is KOI-428. 



Table 3. SOPHIE measurements of KOI-428. 



BJD 


RV 


±lcr,.,. 


BIS 


S/N/pix 


-2 455 000) 


[kms-'] 


[kms-'l 


[kms-'] 


@550nm 


392.5085 


6.320 


0.036 


0.064 


15.7 


395.5188 


6.755 


0.069 


0.033 


10.2 


396.4422 


6.627 


0.029 


0.033 


20.2 


398.5034 


6.411 


0.027 


0.009 


21.0 


400.5389° 


6.549 


0.138 


0.032 


12.0 


401.4932° 


6.850 


0.166 


0.042 


10.9 


425.4872 


6.489 


0.049 


0.073 


12.0 


431.3489° 


6.695 


0.086 


-0.025 


18.3 


441.3261° 


6.582 


0.084 


-0.009 


15.8 


442.4435 


6.634 


0.034 


-0.025 


15.5 


453.4184 


6.417 


0.053 


0.100 


16.1 



Measurement affected by scattered Moon light. 
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